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Abstract. - Based on the mechanics of the Euler equation at short time, we show that a Recent 
Fluid Deformation (RFD) closure for the vorticity field, neglecting the early stage of advection 
of fluid particles, allows to build a 3D incompressible velocity field that shares many properties 
with empirical turbulence, such as the teardrop shape of the R-Q plane. Unfortunately, non 
gaussianity is weak (i.e. no intermittency) and vorticity gets preferentially aligned with the wrong 
eigenvector of the deformation. We then show that slightly modifying the former vectorial field 
in order to impose the long range correlated nature of turbulence allows to reproduce the main 
properties of stationary flows. Doing so, we end up with a realistic incompressible, skewed and 
intermittent velocity field that reproduces the main characteristics of 3D turbulence in the inertial 
range, including correct vorticity alignment properties. 
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Introduction. — Fully developed turbulent flows are 
omnipresent in Nature (e.g. meteorology) and engineer- 
ing (e.g. combustion). Despite an apparent complexity, 
ijt turns out that these flows exhibit universal statistical 
properties such as the Kolmogorov fc~^/'^-law observed on 
t|he power spectrum [l] and the intermittent nature of lon- 
gitudinal and transverse velocity fluctuations [5] . Modern 
developments of experimental and numerical facilities [3]|3] 
have furthermore underlined the peculiar and universal ge- 
ometry of turbulence. 

' From a mathematical viewpoint, our knowledge of 
Navier-Stokes and Euler equations is still poor and this 
rtiakes the situation far from clear if we look for a rele- 
vant model of turbulent flows. Indeed, one can prove the 
existence of dissipative solutions of Euler equations and 
the associated inertial dissipation process has been stud- 
ied [5]. As a consequence the introduction of viscosity is 
not the only way we have to build realistic model of tur- 
bulent flows. Some progresses have been achieved recently 
in the understanding of the universality of the small scales 
of turbulence by studying the Lagrangian dynamics of the 



velocity gradient tensor Ai 



djUi 



where u denotes the 



velocity (see [B] and references therein). While this ap- 
proach takes completely into account the local interactions 
governing the dynamics of A, it requires closures for both 



the pressure Hessian and viscous term. As far as we know, 
provided closures miss at least partially the nonlocal na- 
ture of pressure [B]. An alternative approach would be 
devoted to consider the spatial distribution of the vecto- 
rial velocity field. Very few theoretical works have focused 
on this difficult, although important, aspect. Nonetheless, 
it has been shown that taking into account numerically 
the short time advection of fluid particles for all scales of 
motion allows to build a realistic velocity field that shares 
many properties with empirical turbulence [7]. Unfortu- 
nately, incomprcssibility has to be imposed at every scale, 
leading to a complicated construction method of this field, 
making it not fully explicit. Another route would be de- 
voted to directly proposing a wide class of intermittent 
vectorial fields [5] . But the main limitation of such mod- 
els is that the distribution law of the velocity increments 
is symmetric, so that (following the 4/5 Kolmogorov law) 
such fields do not exhibit energy dissipation. As it is dis- 
cussed in |Hj , it appears that it is a very intricate issue to 
construct an incompressible random field with dissymmet- 
ric increments and non zero dissipation. This issue is the 
subject of the proposed article and while in the symmet- 
ric case explicit analytical calculations can be performed, 
they are out of reach in the present case and consequently 
we make use of numerical simulations to study the model. 
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In this letter, based on former works [HIIH], we propose 
a stochastic method to build an incompressible, skewed 
and intermittent velocity field. This method is motivated 
by the early stage mechanics of the Euler equation dur- 
ing which vorticity is stretched by the local deformation, 
whereas early advection by the large scale velocity is ne- 
glected. We will sec that such a Recent Fluid Deformation 
(RFD) closure [B] leads to an incompressible differcntiable 
velocity field which reproduces well known facts of empiri- 
cal turbulence, namely the teardrop shape of the RQ plane 
and a skewed probability density function (PDF) for the 
longitudinal gradients. Unfortunately, it is shown numer- 
ically that this field is not skewed in the inertial range, 
leading to vanishing mean energy transfer through scales, 
and furthermore, alignment properties of vorticity devi- 
ate from empirical findings. We then show that a slight 
modification of this field, inspired by the multifractal phe- 
nomenology, is able to reproduce the main properties of 
stationary turbulence. 

Recent fluid deformation closure. Let us now 

introduce a flavor of Euler dynamics in the picture. The 
Euler equation writes: 



It + (u.v)u 



dt 

V.u = 



-Vp 



(1) 



It is classical to introduce the vorticity field uj{t, x) = V A 
u(i,x), and take the curl of ([T|) to eliminate the pressure 
and get the Beltrami equation: 



dt ^ ' 



(u.V)u; , 



(2) 



which together with the system V A u = a; and V.u = 
gives a closed equation in a;(t,x). If vorticity vanishes at 
infinity, the solution of this system is given by the classical 
Biot-Savart formula: 



u(t,x) 



1 

An 



Auj{t,y)dy 



(3) 



In what follows, we shall suppose that we have a smooth 
solution of the system ([2]), ([3]) with initial data Wq. Then, 
it will be convenient to introduce the associated La- 
grangian flow X(t, x) defined by the ordinary differential 
equation "^^^^'^^ = u(i, X(i,x)) and X(0,x) = x. Using 

X(i, x), it is easy to see that ([2]) then writes 
(Vu)a;(i, X(t, x)) or equivalently: 



X. 

du;(f,x(t,x)) 

dt 



da;(t,X(t,x)) 
Jt 



= Sa;(t,X(t,x)) 



(4) 



where S is the deformation rate tensor defined by the split- 
ting of the tensor Vu into antisymmetric and symmetric 
parts Vu = io) A . + S. Let us now focus on the short time 
evolution of the system (U) . Since we suppose that the so- 
lution is regular, we can linearize Q in the neighborhood 
of zero, replacing thus S by So (the strain associated to 
the initial vorticity <jJq), which gives: 



a;(t,x) «e*^°a;o(x-tuo(x)) 



(5) 



using the fact that X(t, x) w x-l-tuo(x). In a first step, we 
will neglect the advection of the vorticity by the velocity 
field and only consider the stretching of the vorticity by 
the initial strain tensor So, which gives at time t: 

u(t,x) = -^J Ae*«o('')a;o(y)dy . (6) 

Starting with the Biot-Savart formula, classical calcula- 
tions [9j[l0] give: 



So(y) 



-P.V. 



{y - (t) (g) [{y ~ a) A LOaia-)] 



+ 



Stt 



[(y-cr)Aa;o(cr)]®(y 



dcr 



(7) 



where the integral is understood as a Cauchy Principal 
Value (P.V.) and the tensor product, i.e. x(8)y = Xii/j. 
Now, it is tempting to introduce in formula ([5]) a random 
field Wo which is divergence-free, homogeneous, isotropic, 
Gaussian and with K41 scaling, that is formally [B|: 



a;o(x) = 



where dW(y) = idWi(y),dW2iy),dW3{y)) is the stan- 
dard vector white noise on R^. A more straightforward 
way to do this is to take for ujq (y) the white noise d'W{y) 
and only change to appropriate values the exponents of 
the denominators in the kernels giving u(x) (|x — y]^"^ 
is replaced by |x — y|^(2 + 3)) and So(x) (|x — y|~^ by 
|x - y|~^, with /3 = 2 + I -I- |). Notice that now, the inte- 
gral in the modified ([7]) is no more a principal value. These 
considerations lead finally to define the random field: 



-(PL 



(x - y) A e*^«(y)dW(y) , 



(8) 



with 



(y - t) ® [(y - cr) A dW{(T)] 



1/9 



ly-*^!: 

[(y -a) A dW{(T)] ® (y - (t) 



where, in order to get mathematically well deflned in- 
tegrals, we have introduced both a large scale cutoff 
(/3i(x — y) in the definition of u and So, and a small scale 
regularization e: |x|e = 6^ * |x| (* stands for the convolu- 
tion product) , where is a radially symmetrical function, 
such that /jj3 6'(x)dx = 1 and 6',(x) = ^6'(x/e) 0. 

Numerical study. — We would like now to study the 
statistical properties of the velocity field defined by Eq. 
([5]). Analytical formulas are difficult to obtain, thus we 
will focus on numerical simulations. To do so, one has to 
choose the short time scale t = t. It is easy to check that 
the variance of the matrix So (defined as (tr So)) goes to 
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Fig. 1: Numerical simulations of the process given in ((S)). (a) 
PDF of longitudinal (solid line) and transverse (dashed line) 
velocity gradients for the N = 512 case, (b) Skewness (S) 
and flatness (F) of longitudinal (open symbols) and transverse 
(filled symbols) for the three resolutions: A'' = 128 (o), A'^ — 
256 (□) and — 512 (o). (c) Contour plots of the logarithm of 
the joint probability of the two invariants of A (A'^ = 512 case) 
non-dimensionalized by the average strain Q* = Q / {SijSij) 
and R* = R/ (SijSij)^^^ . The thick line corresponds to the 
zero discriminant (Vieillefosse) line. Contour lines correspond 
to probabilities 10"", 10"^, 10"^, 10"\ 1. (d) PDF of the cosine 
of the angle between vorticity and the eigenvectors of the 
strain (see text) associated to three eigenvalues Ai (dashed- 
dot), A2 (solid) and A3 (dashed). 



infinity as the small scale parameter e goes to zero. So we 
take for r the local normalizing value r = (tr Sq)~^/^, in 
the spirit of the RFD closures provided in Ref. [Hj . 

The simulation is performed in a 1-periodic box 
with collocation points. The infinitesimal vol- 

ume is given by dV — dx^, with dx = 1/N. We 
choose as a regularizing function and large-scale cut- 
off the isotropic normalized Gaussian function V3l(x) = 

6'l(x) {^y^^ cxp{-6\x\^ / L^). This allows to com- 
pute analytically the regularized norm of a vector x, 

-6 I"!'' 

useful for numerical purposes: IxL = — l=e + 

^|x| -I- ilj^^ erf ^1^1 ^ ^ with erf the error function. The 

small-scale cut-off is chosen as e = 2dx and the large one 
as L = 1/2. The kernels of the form x/|x|g entering in 
Eq. ([HI), with a=|-|-|ora = /3, are estimated in the 
physical space in a periodic fashion. White noise compo- 
nents dWi, of zero mean and of variance dV, are generated 
in the physical space using a standard random Gaussian 
generator. Convolution products are then performed in 
the Fourier space. The matrix exponential is evaluated at 



each point of space using a Fade approximant with scaling 
and squaring [TT]. We choose N = 128,256,512. Results 
are displayed in Fig. [TJ 

In Fig. [TJa), we represent the longitudinal and trans- 
verse velocity gradient PDFs for the N = 512 case. We 
see indeed that the longitudinal PDF is skewed, but not 
the transverse one (for symmetry reasons). To further 
characterize the structure in scale of this velocity field, 
we represent in Fig. [TJb), the dependence on the scale 
e of the Skewness S = {{5tuf)/ {{5tuff^ and Flat- 
ness F — {{S^u)'^) / {{SiuY)'^ of the velocity increments 
5iU = u{x + £) — u{x), in both the longitudinal (open 
symbols) and transverse (filled symbols) cases. We see 
that S vanishes and F is consistent with a Gaussian pro- 
cess (i.e. i<" = 3) in the inertial range. This means 
that the weak non-Gaussianity observed on the veloc- 
ity gradients does not survive in the inertial range for 
the velocity increments. To further characterize the local 
structure of this field, we represent the joint probability 
of two important invariants of the velocity gradient ten- 
sor, namely Q = -itr(A^) and R = -itr(A^). This 
so-called RQ-planc has been extensively studied experi- 
mentally and numerically (see [31I3] for comparisons). As 
in empirical data, the RQ-plane is elongated along the 
right-tail of the Vieillefosse line, showing predominance 
of both enstrophy-enstrophy production (upper-left quad- 
rant) and dissipation-dissipation production (lower-right) 
regions. Finally, an important nontrivial property of 3D 
turbulence is the preferential alignments of vorticity with 
the intermediate eigenvector of the deformation [31I3] . We 
represent in[IJd), the probability density of the cosine of 
the angle between vorticity and the eigenvectors e^; of the 
deformation 6 = {u),ex.) with Ai < A2 < A3. We first see 
that the vorticity is preferentially orthogonal to e^^ , as in 
empirical data [5] . It has been observed that the vorticity 
gets preferentially aligned with ex^ , as modeled in [B] . We 
see in Fig. [ijd) that the opposite is observed in our syn- 
thetic field, namely, vorticity gets preferentially aligned 
with . In the following, we will sec that including mul- 
tifractality will allow us to predict, among other features, 
correct alignments. 

Including multifractality. In the first part, we 
have deformed at short times a K41 incompressible Gaus- 
sian field by the deformation part of the Euler flow. As we 
have seen, such a velocity field (Eq. (O) is not intermit- 
tent and moreover, velocity fiuctuations are not skewed in 
the inertial range, i.e. there is no mean energy transfer 
across scales. A first idea would be to iterate the con- 
struction of this field several times and look for a fixed 
point if it exists. Preliminary simulations indicate that 
iterating this construction makes intermittency grow, al- 
though it is not clear if the obtained field is scale invariant 
and still, vorticity alignments are not correct (data not 
shown). An interesting development would be to apply 
this construction at each scale of the flow, in the spirit 
of Ref. [7|. This remains to be explored. Anyway, at 
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this stage, we are missing a basic property of turbulence, 
namely the existence of long range correlations. Indeed, 
the Russian school showed that the dissipation field is cor- 
related over the large integral length scale of the flow, as 
it has been found experimentally (see |12j and references 
therein). Similar observations have been made on the ac- 
celeration in Lagrangian turbulence |13H15| : acceleration 
is correlated over the Kolmogorov time scale, whereas its 
magnitude is correlated over the integral time scale. In 
particular, it has been proposed in jl3j a multifractal ran- 
dom walk able to reproduce this very peculiar property. 
We propose in this part a generalization of this ID inter- 
mittent model of turbulence consistent with the vectorial 
field structure of Eulerian turbulence. 

Following previous works [T6H2T] . it has been proposed 
in general ideas leading to intermittent vectorial fields 
that reproduce the correlation structure of the dissipation 
field. The main underlying idea behind building up an 
intermittent (scalar or vectorial) field is to take the expo- 
nential of a noise correlated logarithmically over a large 
scale L. For turbulence, a way to achieve this is to con- 
sider incompressible homogeneous, isotropic and intermit- 
tent velocity fields of the form 



1 



u,(x) = -- / ^L(x-y) ^ y , Ae^-(y)-^'dW(y) 
47r7 |x-y|| + ^ 

(9) 

where the scalar field X^{y) is defined by the following 
scalar product 



^.(y) = A / 



y cr 



\y-(T\<L |y - cr|, 



— .dW'(cr) 

+1 ^ ^ 



(10) 



Here, the vectorial white noises dW and dW are inde- 
pendent. The constant Ce is chosen such that the ve- 
locity field converges to a non trivial field u as e goes 
to zero [5]. It can be shown rigorously that this vectorial 
field is intermittent with structure function exponents, i.e. 
M,^(|r|) = (|(u(x + r)-u(x)).^|9) ^ |r|C in the longitu- 
dinal case and Mj(|r|) = (|(u(x-t-r)-u(x))A|^|«) - |r|C 
in the transverse one, behaving as a quadratic function of 
the order gr, namely = = 9/3 - 2tt\^ q{q ~ 2). The 
free parameter A entering in the definition of the scalar 
field (Eq. (|10p ) is called the intermittency coefficient. 
It can be easily shown in the limit e ^ |y| — > that 
(X,(y)X,(0)) ~ 47rA2ln(i/|y|). We are indeed taking in 
the model given by the formula Eqs. (|9]) and ((T0| the 
exponential a Gaussian noise correlated over the integral 
length scale L in a logarithmic fashion. Unfortunately, at 
this stage, for symmetry reasons, the velocity increments 
do not exhibit asymmetry and thus, their skewness van- 
ishes for any scale. However, it is proposed in [8] another 
compressible version of this velocity field that is shown to 
provide skewness as soon the very same noise c?W enter- 
ing in Eqs. ^ and ([TU]) was used in the construction. It 
still remains to build up an incompressible and skewed in- 
termittent velocity field that reproduces furthermore non 



trivial geometrical properties such as the asymmetry of the 
RQ plane and the alignments of vorticity with the eigen- 
frame of the dissipation. To do that, one has to consider, 
as it is suggested by the early time advection of fluid parti- 
cles (Eq. ([5])), the exponentiation of a tensor field instead 
of a simple scalar field. 

Motivated by the RFD approach of the first part of this 
article and the explicit intermittent velocity field shown in 
Eqs. © and ([TU)) . we choose to modify directly the field 
given in in order to get a multifractal velocity field. To 
do so, one needs to introduce this unknown intermittency 
parameter A and to change the exponent /3 entering in 
the associated strain So of ([5]) to /3 = 3/2 -I- 2 = 7/2 in 
order to impose logarithmic long range correlations over 
the integral length scale L. Accordingly, we consider the 
incompressible field: 

S,(x) I yPL(x-y) """T^. Ae^(y)dW(y) (11) 

where S is a tensorial Gaussian log-correlated noise, in- 
spired by ([S]) and (|TU)) , of the form 




\y-CT\<L 



(y - t) ® [(y - cr) A rfW((T)] 



,7/2 



[(y - cr) A dW(cr)] ® [y - a) 
U,_^|V2 



(12) 



It can be shown that, for instance, the di- 
agonal components correlation is of the form 
(5'ii(y)5'ii(0)) ~ |A2ln(L/|y|) and, for off-diagonal 
components, (5'i2(y)S'i2(0)) ^ 2A^ln(L/|y|) in the limit 
£ ^ |y| ^ 0- In the tensorial case, calculations are 
difficult and the velocity field (|11[) is expected to be 
asymptotically multifractal [8] with a quadratic structure 
exponent, i.e. for the longitudinal case d^Q^ /dq^ = — cA^. 
A rigorous derivation of the constant c is still missing; we 
will present latterly that numerics show that C<f ~ 1 and 
c « 1 (see the flatness in Fig. [D^b) and the discussion 
of Fig. [3]). The intermittency coefficient A is chosen as 
A^ = 0.025 on empirical grounds [^1^ . 

Numerical results. — We display in Fig. [5] the re- 
sults of simulations of the process given by (|lip with res- 
olutions N = 256, 512, 1024. We first see in Fig. HJa) 
(resp. (b)) the typical continuous shape deformation of the 
longitudinal (resp. transverse) velocity increments PDFs 
characteristic of intermittency and turbulence (see Rcf. 
|H1[53]). The dissymmetry of PDFs in the longitudinal 
case should be noted. In Fig. [H^c), the obtained RQ- 
plane is realistic of a fully developed turbulent fiow. We 
reproduce in Fig. ^d) the skewness and flatness of ve- 
locity increments. First Flatness values are much bigger 
than in the first case and large length scales are popu- 
lated by intermittency. A rough power law is obtained of 
slope —0.1 showing that c « 1 in the longitudinal case 
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Fig. 2: Numerical results of the process given in (|lip. PDFs 
of longitudinal (a) and transverse (b) velocity increments Seu 
(N — 1024), scales £ are logarithmically spaced between dx and 
L. (c) RQ-plane, as in Fig. [Tfc), for A'' = 1024, contour lines 
correspond to probabilities lO'^'^, 10"^ lO"^'^, 10"\ lO"'^, 1. 
Scale dependence of the Skewness (S) and Flatness (F) of the 
velocity increments, as in Fig. [Hb): = 256 (o). A*" = 512 
(□) and iV = 1024 (o). (e) PDF of the cosine of the angle 
between vorticity and the eigenvectors of the strain ex^ (top) 
and (bottom) for the three resolutions. The arrows indicate 
increasing A'', (f) PDF of the cosine of the angle 9 between 
vorticity and the eigenvectors of the strain as in Fig. [2d) for 
N = 1024. 



(see also the following discussion of Fig. O . In the same 
spirit, skewness of longitudinal increments is non-zero at 
any scale in a realistic way. Finally, let us focus on align- 
ment properties of vorticity. In Figs, ^e) and (f) are dis- 
played various PDFs of the cosine of the angle between the 
vorticity and the eigenframe of the deformation. We see 
that as the resolution N increases (Fig. Hl^e)), the vorticity 
gets preferentially aligned with the intermediate eigenvec- 
tor, whereas it gets uncorrelated with the direction of the 
eigenvector associated to the most extensive eigenvalue. 
In Fig. ^f) is reproduced the alignments properties with 
the three eigenvectors of the deformation of the N = 1024 
case, which is realistic of a turbulent flow. We can see 
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Fig. 3: Intermittent characteristics of the process given in (|11|) 
for the A'' = 1024 case, (a) (resp. (b)) Relative behavior of 
the longitudinal (resp. transverse Mj) structure func- 
tions with respect to the corresponding third order moment 
(resp. MJ). Orders considered are represented with dif- 
ferent symbols: o [q — 2), D {q = 4), A (q — 5) and o 
{q — 6). Straight lines correspond to a least-square fit pro- 
cedure. Curves are vertically arbitrary shifted for clarity, (c) 
Scale dependence of the corresponding third order structure 
function M3 /a^ (lower) and /a"^ (upper), with cr^ = 2(uf). 
Dashed line corresponds to a unit slope line, (d) Corresponding 
sets of exponents C,q (o) and (^J (□) assuming that the third 
order structure functions behave as a linear function of the 
scale in the inertial range, in both the longitudinal and trans- 
verse cases. The dashed line corresponds to the K41 prediction 
C,q ~ q/3 and solid curves to quadratic fits (see text). 

here that including the multifractal phenomenology into 
the linearized form (Eq. ([8])) allows to predict surprisingly 
both skewness and peculiar alignments of vorticity. A full 
theoretical explanation of these numerical results is still 
missing, some analytics in this direction would help to un- 
derstand these facts. Let us now focus precisely on the 
intermittent nature of this velocity field, in particular on 
a possible difference between longitudinal and transverse 
velocity increments. 

On the one hand, the turbulence model proposed in Eq. 
pT|) depends on a single free parameter A, the intermit- 
tency parameter. On the other hand, two types of inter- 
mittency can be identified on a vectorial field, namely the 
longitudinal and transverse ones. From both experimental 
and numerical flows, these two types of intermittency seem 
to be different (see Refs. UniUni and references therein). 
As we mentioned, an exact derivation from the model of 
the corresponding longitudinal and transverse (J sets 
of structure functions exponents is still missing because 
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mainly of the difficulty to deal with the tensorial nature 
of the matrix exponential. In order to precisely quantify 
these sets of exponents, we provide a standard fit proce- 
dure of the behavior in scale of the structure functions, in 
a relative fashion in the spirit of Rcf. [21] . 

In Figs. [2Ia) and (b), we perform a least-square fit of 
the relative high orders longitudinal and transverse struc- 
ture functions. This gives access to both the longitudi- 
nal and transverse sets of exponents, assuming for 
the third order structure a linear behavior with the scale 
£ = \r\. We indeed see in Fig. ^c) that in a first ap- 
proximation, and behave as a linear function of 
the scale. In the sequel, we will thus take ^3 ' « 1, this 
could be exactly imposed if the exponent | + | entering 
in Eq. (jlip is slightly changed. Corresponding structure 
functions sets of exponents are reproduced in Fig. EJd). 
One can see that transverse velocity fiuctuations are more 
intermittent than the longitudinal ones, as it has been 
seen in experiments and numerical simulations [251126] . 
The solid curves of Fig. ^d) are quadratic functions of 
q: Qi.T ^ (1 + 3^L,T^2-) q „ c^'^A^^ . Following the fit 
procedure formerly described, numerics show that w 1 
and c'" « 1.4, showing that transverse velocity increments 
are more intermittent than longitudinal ones. These re- 
sults are consistent with the ones obtained from closures 
of the dynamics of the velocity gradients [57]. 

Conclusion and final remarks. As a conclusion, 
based on prior works [Sl[5], we have built an incompress- 
ible skewed intermittent velocity field that reproduces the 
main characteristics of 3D fully-developed turbulence in 
the inertial range. This includes several non trivial geo- 
metrical properties (RQ-plane and the preferential align- 
ments of vorticity), a non- vanishing skewness of the lon- 
gitudinal velocity increments and a realistic intermittent 
picture of longitudinal and transverse velocity fluctua- 
tions. The model provided in Eqs. and can 
thus be seen as a stochastic representation of the fields 
described by Kolmogorov and Obouhkov [2Sl[in]- To do 
so, we included the multifractal phenomenology to the Eu- 
ler mechanics at short-time. Several remarks can be made 
at this stage. First, this theory contains a free parame- 
ter A chosen to be consistent with experimental findings. 
It would be very interesting to find constraints able to 
lead to a direct determination of A. In this work we make 
use of some well known facts, like RFD closure and inter- 
mittency, to build a rather realistic stationary model of 
turbulent velocity field. But we make no use of the very 
dynamics of the system. We may imagine that such a field 
is close in some sense to an invariant measure for Euler or 
Navier -Stokes dynamics and that this would help us to 
fix the value of the unknown parameter A. Hence, time 
correlations and energy (or enstrophy) budgets are not, at 
this stage, predicted. Finally, this model could also help 
in the understanding of the physics of the pressure Hessian 
p[5n]. We leave these aspects for future investigations. 
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